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Abstract 

We have determined the lattice dynamics of MgSiOs perovskite (bridgmanite) by a combination 
of single-crystal inelastic x-ray scattering and ab initio calculations. We observe a remarkable 
agreement between experiment and theory, and provide accurate results for phonon dispersion 
relations, phonon density of states and the full elasticity tensor. The present work constitutes an 
important milestone to extend this kind of combined studies to extreme conditions of pressure and 
temperature, directly relevant for the physics and the chemistry of Earth’s lower mantle. 
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I. INTRODUCTION 


The Earth’s lower mantle, spanning from 670 to 2900 km in depth, represents roughly 
60% of Earth’s total volume. Numerous studies indicate that its chemical composition is 
likely dominated by a perovskite-structured mineral phase, with MgSiOs being its main 
constituent, with Fe and A1 as minor elements^^^. Accordingly, magnesium perovskite is 
likely the most abundant mineral on the Earth. However, irrespectively of its abundance, 
the natural occurrence of MgSiOs perovskite was reported only very recently^, in the Tenham 
L6 chondrite, a shocked meteorite, as the mineral is metastable under ambient condition 
and unpreservable during the exhumation of terrestrial mantle rocks. After the discovery, 
MgSiOs has been named bridgmanite^. 

The determination of the physical properties of bridgmanite is of fundamental geophysical 
interest. Indeed, while the analysis of seismic records allows the determination of compres- 
sional and shear velocities in situ, the establishment of a compositional and mineralogical 
model for the mantle requires the knowledge of the mineral’s equation of state and its elastic 
and thermodynamic properties at the relevant pressure (P) and temperature (T) conditions. 
Structural studies at high pressure and high temperature by X-ray diffraction on samples 
compressed in a laser-heated diamond anvil cell have been extensively performed^^^ and 
provide important insight into the pressure-volume (P-V) and pressure-volume-temperature 
(P-V-T) relationships. The elastic moduli of MgSiOs perovskite at ambient conditions have 
been measured on both single-crystal and polycrystalline samples®’®. The aggregate shear 
modulus has been measured up to 96 GPa^®. Numerous ab initio calculations have been 
performed as well, with the aim not only to validate the experimental structural studies, 
but also to predict thermodynamic and elastic properties, phase stability, and effects of 
anharmonicity via the calculation of the phonon dispersion^^“^^. Noteworthy, despite this 
large body of theoretical calculations, there is not a single experimental phonon dispersion 
study. This is largely explained by the fact that single crystals of sufficiently good quality 
are only available in very small quantities. This excludes the utilization of inelastic neutron 
scattering measurements, but it presents a unique opportunity for inelastic x-ray scattering 
(IXS) with meV energy resolution. In fact, IXS offers the possibility to study samples in 
diamond anvil cells (DACs), thus enabling combined pressure-temperature studies at Mbar 
pressure and temperatures up to 1100 


3 


In the present work we utilize IXS to determine the phonon dispersion of single crystal 
MgSiOs along the main symmetry directions at ambient conditions. The experimental re¬ 
sults are confronted with state-of-the-art ab initio calculations. We benchmark our results 
by comparing the derived elastic moduli with Brillouin light scattering measurements and 
the energies of Raman active modes with experimental values available in literature. 


II. EXPERIMENTAL DETAILS 

The bridgmanite crystals were synthesized starting from 85 mol% MgSiOs -|- 15 mol% 
Mg 2 Si 04 mixture as a silicate source. H 2 O was used as a solvent, and the bulk water content 
of the sample was about 13 wt%. Synthesis conditions were P = 24 GPa and T = 1500 °C, 
with heating duration of 20 min using a Kawai type 5000 ton multianvil press at Okayama 
University. Details of the synthesis procedure of bridgmanite can be found elsewhere^^. The 
water content is estimated to be about 140(52) wt ppm using polarized IR spectra of thin 
sections perpendicular to a, b, and c axes of bridgemanite. 

Optically transparent samples were characterized by single-crystal X-ray diffraction using 
an Xcalibur diffractometer (Rigaku-Oxford Diffraction) using Mo Ka radiation (MoKai 
0.70936 A, MoKq !2 0.71359 A) and a CCD detector. Most of the samples proved to be 
twins, except one. Diffraction data were collected on this sample (50 x 50 x 100 fim 
in size) via omega scans (66.57 mm detector-to-sample distance, 25 seconds per frame, 1° 
frame width) up to 37° in Bragg theta angle. The crystal structure of MgSiOa bridgmanite, 
orthorhombic, Pnma, consists of a framework of corner-linked Si06 octahedra with Mg 
in an irregular 10-fold coordination. We determined lattice constants a = 4.9436(3) A, b 
= 6.9149(8) A and c = 4.7902(6) A. Data collected were merged obtaining a discrepancy 
factor among symmetry related reflections (internal discrepancy factor Rjnt), Rmt = 3.1%, 
and typical rocking curve widths were of about 0.4°. The structural rehnement was carried 
out using the least-squares program SHELXL^^ in the space group Pnma starting from the 
coordinate of Yagi et Anisotropic refinement cycles converged to a discrepancy factor 
Ri=2.5% for 535 Fo>4ctFo and a goodness of £t of 1.106. Bonds and angles are similar 
to literature data^^. The IXS experiment was carried out at the Inelastic X-ray Scattering 
Beamline II (ID28) at the European Synchrotron Radiation Facility in Crenoble (France), 
utilizing the Si (9 9 9) conhguration, which provides an overall instrument energy resolution 
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of 3 meV at 17 794 keV. X-ray beam was focused to a spot size at the sample position of 30 
(horizontal) by 80 (vertical) /rm^ full-width-half-maximum (FWHM). Direction and size of 
the momentum transfer were selected by an appropriate choice of the scattering angle and 
the sample orientation in the horizontal scattering plane. The momentum resolution was 
typically set to 0.28 nm“^ and 0.84 nm“^ in the horizontal and vertical plane, respectively. 
Further details of the experimental set-up can be found elsewhere^^. 


III. LATTICE DYNAMICS CALCULATION 

First principle calculations were performed within the generalized gradient approxima¬ 
tion to the density functional theory in Perdew-Burke-Ernzerhof parameterization as im¬ 
plemented in the CASTEP code^^. The pseudopotentials were of the norm-conserving type 
and created using the optimized method of Rappe et ah with single projectors'^. Three 
reference orbitals were treated as valence states for magnesium and silicon, two for oxy¬ 
gen. The wave functions are expanded in plane waves up to a kinetic cut-off energy of 
800 eV at special k-points on a 4 x 3 x 4 Monkhorst Pack grid resulting in a conver¬ 
gence for internal forces of <10“^ eV/A. The structure optimization was performed using 
the Broyden-Fletcher-Goldfarb-Shannon method by varying cell parameters and internal 
coordinates^®. The calculated lattice parameters at zero temperature agree within 1.1% 
with the experimental values at 77K^^ (see Table SI in the supplementary information). 
The dynamical matrices were computed on a 5 x 3 x 5 Monkhorst Pack grid employing 
density functional perturbation theory and were further Fourier interpolated in the cumulate 
scheme including all image force constants^®”®®. The well-converged internal forces yield a 
maximum error in phonon energies of <0.1 meV. Sum rule corrections for acoustic phonons 
and Born charges were applied with a maximal correction of 3 meV at F. IXS intensi¬ 
ties were calculated in hrst order approximation assuming the validity of both harmonic 
and adiabatic approximation^'. Thermodynamic properties are evaluated within the quasi¬ 
harmonic approximation. The elasticity tensor was calculated using a least-square £t as 
detailed below. A sphere in reciprocal space containing 485 q-points close to F with |q|< 
0.02 relative lattice units (r.l.u.) was used for this £t and the dynamical matrices on these 
q-points were computed using Fourier interpolation. According to the definition of acoustic 
phonons, the deformation energies can be calculated at each q-point, using the dynamical 
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matrix and the same displacement modnlus and polarization for all atomic sites. A given 
choice of displacement and q-vector determine the deformation tensor. In elasticity theory, 
the deformation energy can be expressed as a linear combination of the elastic constants for 
the given deformation tensor. We impose via a least sqnare fit that the two ways of calcn- 
lating the deformation energy coincide for all the considered q-points and polarizations, and 
obtain in this way the elasticity tensor. 


IV. RESULTS 

Figure 1 illustrates representative IXS spectra at selected reduced momentum transfers as 
specihed in the legends and hgure caption. The experimental spectra are characterized by an 
elastic line, centered at zero energy transfer, and several inelastic features in the energy range 
between 20 and 110 meV, corresponding to acoustic and optical phonons. A logarithmic 
intensity-scale is chosen for better visualization of the high-energy phonon modes, which 
have a much weaker intensity. Solid blue lines in the individual panels represent the results 
of the theoretical calculations of the inelastic contribution. The computed energy transfer 
scale was adjusted by a single factor of 1.05, and, for each q-value, the intensities were scaled 
by the same factor to obtain best agreement with the experimental results. The contribution 
at zero energy transfer, due to elastic and quasi-elastic diffuse scattering, is observed in the 
experiments but is not accounted for by the calculation. 

Figure 2 compares experimental and theoretical IXS intensity maps along three selected 
principal directions in reciprocal space. These give a more global overview on the complex 
phonon dispersion and the intensity evolution as a function of momentum transfer Q. Besides 
a strong contribution from the elastic line close to the (0 4 0) point and some artifacts due to 
the interpolation routine used for the reconstruction of the IXS maps from individual £xed-Q 
IXS experimental spectra, we note a very good agreement in terms of phonon energies and 
intensities. Together with the results displayed in Figure 1 this proves that our calculations 
not only correctly predict the phonon energies (after a modest overall scaling), but also 
the phonon eigenvectors, which govern the intensities^^. The slight underestimation of the 
calculated energies can be mainly attributed to limitations in accuracy of the exchange 
correlation functional within the general gradient approximation (see Ref.^® for a detailed 
discussion). 
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TABLE I. Phonon energies of the Raman active modes MgSiOs in meV. The 
values are scaled by 1.05 and compared to experimental values at 290K^^. 
classified by very strong (vs), strong (s), weak (w) and very weak (vw). 


calculation 

experiment 

29.9 

s 

31.0 

s 

30.9 

s 

31.6 

s 

34.9 

s 

35.0 

s 

35.9 

vw 

35.7 

vw 

36.2 

vw 

- 


41.8 

s 

41.8 

s 

42.5 

w 

42.4 

w 

46.2 

s 

45.8 

s 

47.0 

vs 

47.4 

vs 

49.2 

s 

48.6 

s 

54.2 

w 

- 


54.5 

vw 

- 


62.1 

s 

62.9 

s 

62.4 

vs 

62.1 

vs 

63.5 

s 

- 


63.7 

s 

67.7 w 


66.6 

vs 

67.1 

s 

77.0 

w 

- 


78.2 

w 

- 


81.7 

s 

- 


82.6 

s 

82.6 

s 

99.3 

vs 

- 


103.7 

vw 

- 


106.4 

w 

- 



ab initio calculated 
The intensities are 
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FIG. 1. Experimental (open symbols with corresponding error bars) and calculated (blue lines) 
IXS spectra of single crystal MgSiOa at selected reciprocal space points (a) X, b) Y, c) Z d) U 
and e) R. The ab initio calculated intensities were convoluted with the experimental resolution 
function and the calculated energy transfer scaled by 1.05. 

In Figure 3 we display the MgSiOs phonon dispersion along selected high symmetry di¬ 
rections. For the low-energy part we superimpose the experimentally determined phonon 
energies. These were extracted by htting a set of Lorentzian functions convoluted with 
the experimental resolution function to the IXS spectra, utilizing a standard minimiza- 
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FIG. 2. Experimental (left panels) and calculated (right panels) IXS intensity maps along the 
indicated directions on a logarithmic intensity scale. The experimental maps consist each of eleven 
spectra linearly interpolated to 200 Q-points. The calculated intensities are convoluted with the 
experimental resolution function and the energy transfer is scaled by 1.05. 
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FIG. 3. Calculated phonon dispersion relations (solid lines) along selected high symmetry direc¬ 
tions. The left panel is a zoom of the low energy part; the right panel covers the full energy range. 
Experimentally determined phonon energies are plotted as blue dots, the calculated energies are 
scaled by 1.05. 

tion routine. The agreement with the calculation is excellent. Inspection of the complete 
phonon dispersion - which extends up to 120 eV - reveals that the correct assignment of 
an experimentally observed phonon mode becomes increasingly more difficult due to the 
high density of phonon branches. Nonetheless, the already mentioned remarkable agree¬ 
ment between experiment and theory (after a mere rescaling) makes us confident that the 
calculation correctly describes the lattice dynamics of MgSiOs even at the higher energies 
not experimentally covered. The phonon dispersion of the optical branches is rather complex 
and shows pronounced splitting of longitudinal optical and transvers optical branches. The 
splitting is due to large Born effective charges. MgSiOs exhibits 24 Raman active modes 
of which 14 have been determined by experiments^^. The calculated energies are given and 
compared to experiments in Table I. The values were assigned by careful comparison of cal¬ 
culated and measured intensities. The mode at 99.3 meV should be visible in experiment, 
but is beyond the measured energy transfers of reported measurements. The agreement 
between experiment and calculation is very good with the only exception of the mode with 
a calculated energy of 63.7 meV. 
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FIG. 4. Calculated phonon density of states of MgSiOa. The partial density of states are given per 
species. The energy is scaled by 1.05. 

The so-validated theoretical results were then used to derive the partial and total phonon 
density-of-states (DOS), shown in Figure 4. The partial contributions per species indicate 
that the vibration of oxygen atoms dominates the hrst two low energy peaks in the DOS. 
Between 67 and 82 meV all modes are almost exclusively due to oxygen. We further note 
a gap between 100 and 102 meV. The total DOS allows the derivation of important ther¬ 
modynamic properties via integral equations^^. Some relevant thermodynamic quantities 
calculated within the quasi-harmonic approximation for room temperatures can be found in 
the supplementary information (see Table S2). 

The full elasticity tensor is given in Table II. The elastic moduli were determined using 
a least square fit to acoustic phonons without the need of symmetry braking displacements. 
For comparison we also employed the traditional method exploiting the stress-strain relation 
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TABLE II. Elastic moduli in GPa derived from the present ab initio calculations, by fitting the 
elasticity tensor to acoustic phonons at low q (DEPT). Experimental values from BLS are reported 
for comparison®. 


Cll C22 C33 C44 C55 C66 C12 C13 C23 

DEPT 511(2) 426(2) 494(2) 176(1) 151(1) 193(1) 149(1) 109(1) 137(1) 

BLS 528 456 481 182 147 200 146 125 139 


from 18 distorted structures, which leads to comparable values. The calculated values are 
in good agreement with literature values determined by Brillouin light scattering (BLS)®. 

The applied method provides a useful tool to derive the elasticity tensor directly from 
the lattice dynamics. It prohts from the acoustic sum rule corrections, which overcome 
numerical issues due to the hnite grid sampling of reciprocal space and hnite number of 
plane waves. As further advantage, in contrast to stress-strain calculations using distorted 
structures, all structural symmetry operations can be exploited to reduce the number of 
k-points for the Brillouin zone sampling. The calculations are thus much easier to converge. 


V. SUMMARY AND PERSPECTIVES 

Our experimental study combined with parameter-free model calculations provides the 
complete lattice dynamics description of MgSi 03 -brigmanite, the main constituent of the 
Earth’s lower mantle and one of the most important Earth’s minerals. More in general, our 
work clearly demonstrates that the phonon related properties of many geophysically relevant 
materials can nowadays be obtained without the need to record the phonon branches one- 
by-one along all the high symmetry directions. We have shown that the determination of a 
representative sub-set of experimental data is sufficient for the reliable determination of the 
full lattice dynamics. Consequently, very important thermodynamic and elastic properties 
can be derived with high accuracy. In particular, the knowledge of the full phonon dispersion 
is crucial for validating any physical quantity related to the lattice dynamics. Furthermore, 
when extended to high pressure and high temperature conditions, it allows the validation of 
the phase stability proposed on the basis of energy minimization. It is worth noting that the 
availability of a single crystal, while beneficial, is not mandatory, and the very same concept 
can be applied to polycrystalline materials as was demonstrated in the case of stishovite^^. 
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The here-reported study at ambient conditions marks the starting point for a comprehensive 
investigation of the stable phases in Earth’s mantle, including Mg 2 Si 04 (forsterite, spinel, 
/3-phase), MgSiOa (orthopyroxene, ilmenite, perovskite), and Si02 (stishovite, seifertite) at 
relevant pressure and temperature conditions of the Earth’s lower mantle. Such a study is 
technically feasible as the x-ray beam at ID28 can be focused down to 17 x 10 /im^ with a 
flux reduction of only 25-30% compared to the optical conhguration used in the present 
experiment. This allows working efficiently with samples of lateral dimensions of few tens 
of microns. Once considering samples down to a thickness of 15 - 20 pm, compatible with 
experiments in a diamond anvil cell (DAC) up to Mbar pressures, the expected loss in 
the IXS signal is 5-7 with respect to the present experiment. This would lead to typical 
acquisition times per spectrum of about 4h, and a complete set of data in a few days. It 
is worth noting that modern IXS spectrometers are equipped with several crystal analyzers 
(nine in the case of ID28); thus IXS spectra at several different momentum transfers are 
recorded simultaneously. 
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SUPPLEMENTARY INFORMATION 


TABLE SI. Lattice constants of MgSiOa perovskite (space groupe Pnma). 
Calculation [A] Experiment (77K) [A]^^ Our experiment (297 K) 


a 

4.979 

4.930 

4.9436(4) 

b 

6.973 

6.896 

6.9149(8) 

c 

4.828 

4.773 

4.7902(6) 


TABLE S2. Calculated thermodynamic properties per unit cell at 298K. Vibrational contribution 
to the total energy E and free energy E; Entropy S and heat capacity C„. 

E (eV) E(eV) S (J/mol/K) C„ (J/mol/K) 


0.476 -0.252 244.4 332.9 
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